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We study two a priori unrelated constructions: the spectrum of edge modes in a band topological 
insulator or superconductor with a physical edge, and the ground state entanglement spectrum 
in an extended system where an edge is simulated by an entanglement bipartition. We prove an 
exact relation between the ground state entanglement spectrum of such a system and the spectrum 
edge modes of the corresponding spectrally flattened Hamiltonian. In particular, we show that 
degeneracies of the entanglement spectrum correspond to gapless edge modes. 
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I. INTRODUCTION 

Topological phases of matter, occurring for example in 
the quantum Hall effect, cannot be distinguished using 
a local order parameter. Although they can sometimes 
be characterized using, for example, topological ground 
state degeneracy or the existence of gapless edge modes, a 
full understanding is still lacking. One worthwhile direc- 
tion is to use the information-theoretic concept of entan- 
glement to characterize topological phases, as was shown 
first in [H, [2| , who related the universal sub-leading term 
in the bipartite entanglement entropy (the "topological 
entanglement entropy" ) to the total quantum dimension 
in a gapped anyonic system. Further progress was made 
by Haldane and Li ^ who related the full spectrum of the 
reduced density matrix (the "entanglement spectrum") 
to the conformal field theory (CFT) edge mode spectrum 
in fractional quantum Hall states. Here, our goal will 
be to study topological phases occurring in free fermion 
systems, i.e. band topological insulators and supercon- 
ductors. In such systems the topological entanglement 
entropy can vanish, so one needs a finer method to dis- 
criminate between the phases. We will adhere to the phi- 
losophy of 3 and focus on the entanglement spectrum, 
which we will show does contain more information about 
the various free fermion topological phases. Specifically, 
we will show that the entanglement spectrum of any band 
insulator or superconductor can be exactly reconstructed 
from the edge mode spectrum of the corresponding spec- 
trally flattened Hamiltonian. 

Given a division of a quantum system into two sub- 
systems A and i?, such that the total Hilbert space is 
the tensor product of the subsystem Hilbert spaces, one 
defines the reduced density matrix pA on A by tracing 
out the degrees of freedom of B from the pure ground 
state density matrix (for fermion Hilbert spaces 

one actually has to introduce graded tensor products, 
but we can safely ignore this subtlety). The von Neu- 
mann entanglement entropy between A and B relative 
to the ground state is then defined hy S = —tv pA log pA. 
This entropy is a measure of the complexity of the ground 
state, as seen, for example, in the logarithmic scaling of 
entanglement entropy at criticality [4!|, as well as in the 
study of one dimensional gapped systems, where finite 
entropy matrix product states (MPS) are introduced to 
approximate ground states Q ■ 



As motivated by Q , we take the more general approach 
of studying the entire spectrum of eigenvalues of the re- 
duced density matrix (i.e. the eigenvalues of the Schmidt 
decomposition of the ground state), in the hope that the 
spectrum contains more information than just the one 
number S that can be constructed from it. We restrict 
our attention to free fermion systems, which, though non- 
interacting, include effective Hamiltonians for topological 
insulators, such as HgTe in 2 dimensions 0,i and BiSb 
in 3 dimensions [lo| , as well as superconductors such 
as SrRu, and systems with broken T symmetry such as 
the integer quantum Hall effect. 

There has been a great deal of work done on entangle- 
ment in free fermion systems. First of all, entanglement 
in a Fermi gas is well understood [3, [ll| . The general- 
ization to arbitrary free fermion systems has been stud- 
ied as well |12| - [l4| . and the problem of computing the 
entanglement spectrum has been reduced to diagonaliz- 
ing a matrix of Green's functions. Here we derive this 
formula in the formalism of free Majorana fermions and 
Gaussian states [isl . IT6j . The advantage of this approach 
is that it unifies the treatment of both band topologi- 
cal insulators and superconductors in that the effective 
Hamiltonian, which can include both hopping and pair- 
ing terms, becomes a general quadratic form in the Ma- 
jorana fermions. We use this formula to relate the en- 
tanglement spectrum to gapless edge modes of a sample 
with boundary. Before we go into more detail, we make 
some comments on related work. 



In 17[ , a numerical approach was pursued in the study 



of the entanglement spectrum of a topological super- 
conductor. Specifically, it was shown numerically that 
one can meaningfully distinguish between the weak and 
strong pairing phases of a 2 dimensional p + ip super- 
conductor by looking at the degeneracies of the entan- 
glement spectrum. In this paper we analytically prove 
a relation which generalizes this result. Furthermore, in 
[Tsl . which appeared shortly after this paper, a general 
relation similar to ours was derived. Working in the set- 
ting of topological insulators, [18] show that gapless edge 
modes imply degeneracies in the entanglement spectrum. 
However, they also show that the converse is not true by 
exhibiting a specific model of a 3 dimensional topologi- 
cal insulator in which a Zeeman field gaps out the edge 
modes but does not effect the degeneracies of the en- 
tanglement spectrum, which are protected by inversion 
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symmetry. 

Let us now describe our derivation in more detail. We 
prove a relation between two a priori unrelated construc- 
tions. The first is the computation of the entanglement 
spectrum of the ground state of a band topological insu- 
lator or superconductor with respect to a partition into a 
large but finite region A (such as a disc in 2 dimensions, 
for example), and its complement B. For the second, we 
need the 'spectral flattening trick' [l^, which we use to 
deform the original Hamiltonian iJ to a new one H' that 
retains the same ground state, but has a flat spectrum; 
this is always possible for a gapped Hamiltonian. The 
deformation keeps the gap open and can be done adia- 
batically, showing that the two Hamiltonians are in the 
same topological phase. We now simulate a physical edge 
at the boundary of A by defining to be the restric- 
tion of the spectrally fiattened Hamiltonian to region A 
- that is, we retain only the couplings within A, and dis- 
card the degrees of freedom in the now decoupled region 
B. We prove that the spectrum of this restricted, spec- 
trally flattened, Hamiltonian H'j^ and the entanglement 
spectrum of the ground state can be reconstructed from 
each other. 

Of particular interest are the low energy modes of i/^, 
which, because the bulk is gapped, are edge modes local- 
ized at the boundary of A. One consequence of our for- 
mula is that zero modes of H'j^ correspond to degeneracies 
in the ent ang lement spectrum. An important point, ob- 
served in [13] , is that the edge modes of H'^ might not be 
the same as those obtained from the restriction Ha of the 
original Hamiltonian to A, without first applying a spec- 
tral flattening transformation. Indeed, as was shown in 
[m , it is sometimes the case that the edge modes of Ha 
are gapped, even though there are degeneracies in the en- 
tanglement spectrum (protected by inversion symmetry, 
for example). The crux of the issue, thus, is that the edge 
mode spectrum can change during the spectral flatten- 
ing transformation. However, if there is a symmetry, such 
as time reversal, that protects the gapless nature of the 
edge modes, then all the Hamiltonians along the spectral 
flattening deformation respect this symmetry, and pos- 
sess gapless edge modes. In this case, therefore, we can 
conclude that there are at least as many corresponding 
degeneracies in the entanglement spectrum, in a sense 
made precise below. As [l8| propose, the entanglement 
spectrum may in fact be a more robust characterization 
of a topological phase than the edge mode spectrum (see 
also fig]). 

Before we move to the general analysis, we illustrate 
our point with a simple example. Consider a one dimen- 
sional superconductor, i.e. a fermionic Hamiltonian with 
both pairing and hopping terms. The Hamiltonian is: 

i7 = i ^ (uaiaj + vaiai+i +va]ai+i\ + h.c. (1) 



forces each site into occupation number or 1, depend- 
ing on the sign of w. At |u| — \v\ there is a critical point 
describing a transition into a different phase {\u\ < \v\) 
(in fact, a Jordan- Wigner transform maps this to the 
transverse field Ising model transition). This new phase 
is topological, in the sense that it is characterize d by 
having unpaired Majorana modes at the boundary [20|. 
That is, for chain whose length L is much larger than the 
correlation length ^, there are two ground states, degen- 
erate up to a splitting exp(— L/^). They form a two state 
system obtained from pairing up the boundary Majorana 
modes. However, we will find below that the entangle- 
ment spectrum also contains a signature of the topolog- 
ical phase. Namely, we will see that in the topological 
phase, the multiplicity of all eigenvalues in the Schmidt 
decomposition is doubled (see fig. [T]). 



trivial phase 



non-trivial 
phase 



FIG. 1: Schematic representation of the entanglement spec- 
trum of the Majorana chain. In the nontrivial topological 
phase, there is a two-fold degeneracy in the spectrum. 



II. GENERAL ANALYSIS 

We consider a free fermion system (i.e. band insulator 
or superconductor) in an arbitrary number of dimensions 
d. We picture a tight binding model with short range 
interactions, which could be either hopping or pairing. 
To conveniently work with both, we write each orbital 
in terms of two Majorana fermions, so that N physi- 
cal fermion modes are described with 2A^ operators Cj, 
j = 1, . . . , 2N. The Cj are Hermitian and satisfy the Ma- 
jorana commutation relations {cj,Ck} = 2Sjk. We can 
write N physical fermion creation and annihilation oper- 
ators as 



This is the so-called Majorana chain Hamiltonian [20|, 
and it exhibits two phases. For |u| > \v\, the chemical 
potential dominates, and in the limit of large u/v simply 



an = {C2n + iC2n~l) 
On = 2 - iC2n-l) ■ 



(2) 
(3) 
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The effective Hamiltonian can now be written in terms 
of the Cj'. 

. 2N 

where we interpret j as a collective site and band index. 
The matrix Hji^ is real skew-symmetric, and we assume 
the Hamiltonian ([4]) is gapped. 

We now choose a partition of the system into a sub- 
system A and its complement B. For instance, in two 
dimensions we could define A and B by partitioning the 
lattice into two complementary half-planes. However, for 
convenience we shall assume A to contain finitely many 
lattice sites. Indeed, let A contain m orbitals, described 
by the Majoranas ci, . . . , C2m, and let B be described by 
the remaining Majoranas C2m+i, • • ■ , C2n- 



What is the M matrix ([5]) for the ground state of the 
Hamiltonian (U)? We can determine it from the 2 point 
functions of the ground state, obtained from the formula 




{'ip\cjCk\ii) = -sg'D-{iHjk) ■ (7) 

The sgn function is defined as follows: for a diagonal 
Hermitian matrix Z), it replaces all positive eigenvalues 
with 4-1 and all negative eigenvalues with —1. A general 
Hermitian matrix, such as iH in ([7]), can be diagonal- 
ized with a unitary transformation U: iH = U DU~^, 
and we define sgn (iHjk) = U sgn (Z?) U^^. Thus, after 
comparison to eq. ^ , we see that 

Mjk = -isgn{iHjk). (8) 
B. Exact correspondence 



A. Gaussian States 

Before proceeding, it will be useful to introduce the 
concept of a Gaussian state p^ . Intuitively, a Gaussian 
state is simply the formalization of what it means for a 
(possibly mixed) state to be the ground state of a free 
fermion Hamiltonian. Since free fermion Hamiltonians 
can be diagonalized, one can think of a Gaussian state of 
as a tensor product of independent 2-state system den- 
sity matrices in some orthogonal basis. More formally, 
consider a density matrix p, written as a polynomial of 
the 2N Majoranas Cj in such a way that each cj occurs 
to the power or 1 in each term. The state defined by 
p is said to be Gaussian if, upon replacing the Cj with 
anti-commuting Grassman variables 9j, one obtains an 
expression p that can be put in the form 

P=^exp(^^0^Af0) (5) 

for some real antisymmetric 2^^^ by 2N matrix M . The 
matrix M simply encodes the 2 point correlators of the 
Cj in the state p: 

Mjk ^Tr{picjCk) (6) 

for j ^ k, with Mjj = 0. All higher correlators are 
determined by Wick's theorem. 

The states with which we will be dealing are all Gaus- 
sian. First of all, the ground state of a gapped 
Hamiltonian is Gaussian (to see this, simply bring Hjk to 
canonical block diagonal form). Furthermore, given any 
Gaussian state p of the full system, the reduced density 
matrix pA constructed from it by tracing out the degrees 
of freedom in B is also Gaussian. Indeed, because the 
correlators of pA are the same as those of p, and a state 
is determined uniquely by the set of all its correlators, 
we see that pA is a Gaussian state whose matrix M ^ 
is simply the restriction of that of p. 



Having introduced the necessary formalism, we now 
prove an exact relation between the entanglement spec- 
trum of a general gapped Hamiltonian, and the edge 
mode spectrum of the corresponding spectrally fiattened 
one. Let us first explain the spectral fiattening transfor- 
mation [15] . Given any gapped Hamiltonian ^ we can 
construct a one parameter family of gapped Hamiltoni- 
ans that interpolate between H and 

. 2N 

= i E ^i^^^-^'^' (9) 

where the eigenvalues of i H'^j, are all ±1. The matrices 
that interpolate between Hjk and iZjj. share a common 
eigenbasis; they are all gapped, and leave the ground 
state invariant. Also, as a consequence of the gap being 
open, all of the Hamiltonians in the family are quasi- 
local, with quadratic coupling terms exponentially sup- 
pressed by the distance |15| . 

Because the spectral fiattening transformation leaves 
the ground state invariant, the entanglement spectrum, 
which depends only on the ground state, is the same for 
H and H' . Thus, to prove our relation we will from now 
on just assume H is spectrally fiat, and omit the extra 
superscript in H' . With this assumption, equation (|8]) 
simplifies to 

Mjk = Hjk- (10) 

Equation (fTU|) is the key to proving the correspondence. 
Roughly speaking, we will relate its right hand side to the 
edge mode spectrum, and its left hand side to the en- 
tanglement spectrum. More formally, let Hjk, and Mjk 
denote the restrictions of Hjk and Mjk to A respectively; 
thus Hjk and Mjk are 2m by 2m anti-symmetric matri- 
ces, and they are equal by virtue of fTU]) . The matrix 
Hjk now defines a physical Hamiltonian on region A: 

Ha = ^Y1 ^jkCjCk- (11) 

j,k£A 
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Ha is gapped in the bulk, with spectral gap normalized 
to 1 in our units, but it also has boundary modes, which 
could potentially be gapless. The eigenvalues of iHjk, 
which come in pairs ±Ar, |Ar| < 1, r = l,...,m, re- 
flect these boundary modes, in that the corresponding 
eigenstates are localized near the boundary whenever A,- 
differs substantially from ±1. 

From our discussion of Gaussian states, on the other 
hand, we see that Mj^ is just the M matrix (see eq. 
E]) of PA, the reduced density matrix obtained from the 
ground state IV') ("01 by tracing out the degrees of freedom 
in B. This means that if one defines pA, as in ([5]), to be 
PA with the Majoranas Cj, j = 1, . . . , 2m, replaced with 
anti-commuting Grassman variables 9j , then 

PA = ^ exp (^lo^MO^ - ^ exp (^^^^i^^) ■ (12) 

Since the expression 9^ HO is S0(2m) invariant, we can 
rotate to a more convenient basis. In particular, let H'^ 
be the canonical block diagonal form of H, consisting of 
m blocks of the form 

(-°a/o')' 

Let M be the SO {2m) matrix that block diagonalizes 
H: H'= = MHM-^, and define 9'^ = Efe^^jfe^fc' and 
Cj- = 'EkMjkCk- 

In terms of the rotated variables, ((T^ turns into 

When we expand out the last expression in ([HI) , each 9^ 
occurs to power or 1, and hence 

= ft + ^4.-lC2.) ■ (15) 

From the product form (jl5l) of pA, we can immedi- 
ately read off the entanglement spectrum. Indeed, ([15]) 
shows that the density matrix pA decomposes into m in- 
dependent 2 state systems. The r'th one, where r = 
1, . . . , m, described by 1/2 -f i (Ar/2) c'2r^ic'2r, has eigen- 
values l/2±Ar/2. Hence the full entanglement spectrum 
is described by the set of 2™ Schmidt eigenvalues: 




III. CONSEQUENCES 

We have, in ([T6l) . computed the entire set of 2™ 
Schmidt eigenvalues comprising the entanglement spec- 
trum of the ground state with respect to regions A and 
B in terms of the eigenvalues itA^. These eigenvalues 
determine the edge mode spectrum of the corresponding 
spectrally flattened Hamiltonian ([TT|) . Note that they 
only differ significantly from ±1 when the correspond- 
ing eigenstate is localized near the boundary of region A. 
Also, if Xr ^ ±1, which corresponds to a bulk mode, then 
1/2 ± Xr/2 are close to and 1; hence the contribution 
of such a mode to the entanglement becomes vanishingly 
small as Ar — > 1 (or Xr —1)- This just means that only 
edge modes contribute significantly to the entanglement. 
Also, we see that a zero energy edge mode, correspond- 
ing to Xr = 0, is reflected as a nontrivial multiplicity of 
the entanglement spectrum (|16p . Indeed, having k such 
Ar = results in a multiplicity of 2*^ for the Schmidt 
eigenvalues. 

We note the versatility of our Majorana fermion ap- 
proach: our results hold not only for topological insula- 
tors, but for topological superconductors as well. Fur- 
thermore, because the entanglement spectrum depends 
only on the ground state, our approach gives a way to 
diagnose topological order by looking at only the ground 
state wavefunction. Indeed, we have proved that the 
entanglement spectrum can be reconstructed from the 
edge mode spectrum of the corresponding spectrally flat- 
tened Hamiltonian. For a generic topological Hamilto- 
nian with symmetry protected edge modes, we expect 
that the corresponding spectrally flattened Hamiltonian, 
being in the same phase, has the same symmetry pro- 
tected edge modes. Thus the entanglement spectrum 
is as good as the edge mode spectrum at discriminat- 
ing between the topological phases. However, the above 
statement can fail in special cases: for example, inver- 
sion symmetry can sometimes allow the spectrally flat- 
tened Hamiltonian to have a richer gapless edge mode 
structure than that possessed by the original Hamilto- 
nian [l^. In fact, [l^ and [l^ interpret this discrepancy 
as a virtue, and argue that the entanglement spectrum 
is a more robust measure of topological order than the 
edge mode spectrum. 

It would be very interesting to generalize this corre- 
spondence to interacting systems, as was done, for ex- 
ample, by Li and Haldane in [^] for certain fractional 
quantum Hall states (see also |2l|V 
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